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Abstract 

An empirical study of the bivariate probability distribution characterizing 
a full set of two consecutive price increments x and y for a group of stocks 
at time scales ranging from one minute to thirty minutes reveals asymmetric 
structures with respect to the axes y = 0, y = x, x = and y = —x. All 
four asymmetry patterns remarkably resemble a four-blade mill called market 
mill pattern. The four market mill patterns characterize different aspects of 
interdependence between past (push) and future (response) price increments. 
When analyzed in appropriate coordinates, each pattern corresponds to a 
particular nonlinear dependence between the push and the conditional mean 
of response. Qualitative interpretation of each pattern is discussed. The 
market mill pattern is an evidence of complex dependence properties relating 
past and future price increments and resulting in various types of nonlinear 
correlation and predictability. 
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1 Introduction 



The apparent randomness of stock price dynamics explains a dominant role 
of the random walk paradigm in theoretical finance. The evolution of stock 
prices is indeed close to being random, but a careful statistical analysis re- 
veals a number of characteristic dependence patterns such as correlations 
between simultaneous increments of different stocks, statistically significant 
autocorrelations at intraday timescales, volatility clustering, leverage effect, 
etc., that became a part of a standard toolkit for model - building [TJ |21 El El ■ 
Each of these effects corresponds to some sort of probabilistic dependence 
between lagged and/or simultaneous price increments. 

A general way of quantifying the probabilistic dependence providing a 
full description of the set of quantities under consideration is to specify the 
corresponding multidimensional probability distribution, from which one can 
calculate all quantities of interest such as various conditional distributions 
and their moments. Generically such distributions can depend on simulta- 
neous and lagged price increments, traded volumes, etc. 

Our approach to studying the dynamical patterns in the stock price evolu- 
tion OinilZ] is based on a direct analysis of such multidimensional probability 
distributions. A simplest case we concentrate upon in our studies jSl El U\ is 
that of a bivariate distribution describing the interdependence of two price 
increments in two coinciding or non-overlapping time intervals. Even here 
there is quite a few cases to analyze: the increments can be simultaneous and 
belong to different stocks or lagged and belong to the same stock or different 
stocks, the time lag separating the intervals can vary, etc. The interde- 
pendence between lagged price increments for the same stock was termed 
"horizontal" and the one between simultaneous price increments of different 
stocks termed "vertical" in where both horizontal and vertical dependen- 
cies were considered. In the present paper we focus on the horizontal case 
only. 

Despite of the fundamental importance of the issue, a list of papers on 
the multidimensional probability distributions of stock price increments is, to 
the best of our knowledge, not abundant. A theoretical analysis of bivariate 
distribution of returns in two consecutive intervals in the particular case of 
Levy-type marginals was performed by Mandelbrot jHj, where some interest- 
ing geometric features of this distribution both for the case of independent 
and dependent returns were described. As discussed in 0, the bivariate dis- 
tribution in question can be considered as a "fingerprint" reflecting a nature 
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of the pattern embracing the two consecutive returns. 

Recently an analysis of conditional dynamics of volatility exploiting con- 
ditional distributions was described in JU] ■ Let us also mention direct studies 
of the above- described bivariate distributions for returns in consequent time 
intervals that revealed interesting phenomena such as the "compass rose" 
[TTl IT2*1 ITB*] . The conclusions on whether the discovered patterns are re- 
lated to predictability differed: whereas [TT] found no predictability, in [T3] 
an interesting pattern of wavelet-filtered series corresponding to partial pre- 
dictability was found. The first moment of corresponding conditional distri- 
bution for daily time intervals was studied in . Another line of research is 
an explicit analytical reconstruction of the bivariate distribution in question 
using copulas, see e.g. |T51 EQU UH HE] • 

Any dependence between lagged price increments means that certain 
characteristics of price evolution are predictable. The issues of predictability 
and market efficiency are of paramount importance for theoretical finance. 
The analysis of the present paper is relevant in this context because a key 
quantity characterizing predictability, the conditional mean [TH1 I2HI I2H 1221 
123] . is determined by the conditional distribution of the increment in the 
second interval at given increment in the first one. Let us emphasize that 
the profitability of a strategy given the existence of the nonzero conditional 
mean requires, in each and every case, a separate analysis. A particular ex- 
ample of potential profit from high frequency anticorrelations of stock price 
returns was discussed in pQ. The properties of conditional mean are most 
often discussed in the context of verifying the martingale property of the 
financial series [23 which is, in turn, a generally accepted outcome of market 
efficiency. 

The paper is organized as follows. 

In paragraph 2.1 we describe the market data used in our study and 
elaborate on basic points behind the probabilistic methodology used in its 
analysis. 

Main empirical results of the paper are described in paragraph 2.2. We 
begin with discussing the symmetry properties of the bivariate distribution 
of price increments in two consecutive time intervals, termed push and re- 
sponse. The distribution is shown to be approximately invariant with respect 
to rotations at a multiple of tt/2. This gives rise to a number of relations 
between the probabilities of certain two-step processes. We continue with 
discussing four asymmetries of the bivariate distribution with respect to the 
axes y — 0, y = x, x = and y = —x: asymmetry of conditional response, 
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asymmetry with respect to increment's ordering , asymmetry of conditional 
push and asymmetry of mirror orderings correspondingly. Remarkably all of 
these asymmetries are characterized by the same geometric pattern termed 
market mill. The asymmetries are shown to result from the basic geometry 
of the full distribution described at the beginning of the section. 

In paragraph 2.3 we discuss major outcomes from the market mill pat- 
terns. In particular a nonlinear dependence of the conditional mean response 
on the push clearly demonstrates the nonlinear nature of correlations between 
the two variables in the considered bivariate probability distribution. The 
market mill pattern signals the market inefficiency at intraday time scales 
resulting from temporal interdependencies between consecutive price incre- 
ments. 

We conclude by formulating the main results and outline of the future 
work in paragraph 3. 

2 Asymmetry patterns in price dynamics 
2.1 Data and Methodology 

Our study of high frequency dynamics of stock prices are based on data on 
the prices of 100 stocks traded in NYSE and NASDAQ in 2003-2004 sampled 
at 1 minute frequency 2 . 

A method of analyzing dependence patterns in financial market dynamics 
adopted in the present paper, see also El El , is based on a direct exami- 
nation of multivariate probability distributions and conditional distributions 
constructed from them. Below we shall restrict our consideration to the bi- 
variate distribution of stock price increments in two non-overlapping time 
intervals. 

Let us consider two non-overlapping time intervals of length ATi and AT 2 . 
The time ordering is such that the interval AT 2 follows after the interval A2\. 
Generically the end of the first interval and the beginning of the second can be 
separated by some time interval but in most of cases under consideration the 
second interval immediately follows the first one. In this case we are dealing 
with two price increments Sp(ATi) = p(t) — p(t — A7\) and 5p(AT 2 ) = 
p(t + AT 2 ) —p(t). In what follows we shall denote the price increment in the 

2 All holidays, weekends and off-market hours are excluded from the dataset. The list 
of stocks is given in the Appendix. 
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first interval 5p(ATi) by x (push) and that in the second one Sp(AT 2 ) by y 
(response). Below we consider the case of A7\ = AT 2 = AT and study four 
time scales AT = 1 min, 3 min., 6 min. and 30 min. 

The full probabilistic picture of interrelation between the push and re- 
sponse is given by a bivariate probability density V(x,y). This is the most 
comprehensive quantity from which other quantities of practical interest 
could be derived. In particular, one can compute the conditional proba- 
bility distribution V(y\ x) = V(x,y)/V(x) describing the probabilistic shape 
of response at given push. 

Let us stress that in the present paper we consider the probability distri- 
bution characterizing a full set of price increment pairs for a group of stocks 
so that a set of events (pairs of increments) under consideration unifies all 
subsets of events characterizing individual stocks. Studying the properties of 
thus obtained bivariate distribution V(x,y) is the main goal of the present 
paper and the companion paper ||. Bivariate distributions for individual 
stocks are considered in the companion paper [7|. 

Let us now comment on our choice of price increments and not price 
returns as basic variables. The question of a choice of the "right" variable 
is by no means trivial. The full picture may well be the one using both 
increments and returns at different time scales. At intraday time scales 
the market participants use strategies based on evaluating price moves in 
dollars, while at larger time scales (e.g. weeks) they are rather formulated 
and realized in terms of returns. The particular scheme interpolating between 
price increments at small time scales and returns at large ones is described in 
[3]. Let us emphasize that we have checked that the main results described 
in the paper are valid for normalized returns as well. 

2.2 Empirical results: market mill patterns 

Let us now turn to discussing the structure of push - response relationship 
by considering the market data on price increments in two consecutive 1 - 
minute intervals. As discussed in the previous section, the full probabilistic 
picture of the two-step push - response process is described by the probability 
distribution V(x,y). The two - dimensional projection of \og 8 V(x,y) 3 is 

3 In all plots we are using the log 2 n scale for the probability distribution so that the 
equiprobability lines differ by the factor of 2™. The particular values of n are chosen to 
show the distribution geometry in the most clear fashion. The borders between different 
colors correspond to the equiprobability lines. 
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plotted in Fig. 10 4 . This distribution turns out to have many remarkable 
properties that we analyze below and in the companion papers jHllZl- 

To facilitate the analysis let us sketch the shape of the equiprobability 
levels in Fig. 1, in which the xy plane is divided into sectors numbered 
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Figure 1: Sketch of the equiprobability levels of the bivariate distribution 
V(x,y). Straight lines are drawn to guide the eye. 

counterclockwise from I to VIII. 

The first striking property of the bivariate distribution V(x,y) is its ap- 
proximate invariance with respect to rotations at the multiples of ir/2. In 
terms of sample paths composed by the increments ±d and ±£2 the symme- 
try with respect to rotations at multiples of ir/2 leads to a chain of equalities 5 

V(&, C2) = V{-{ 2 , CO = V(-Ci, -C2) = V(( 2 , -Ci) (1) 

Second, let us stress that naively expected commutativity of price increments 
under addition T-'(Ci) C2) = ^{^2, Ci) is violated by the asymmetries of the 
bivariate distribution V(x, y) analyzed in full details in the next paragraph. 

4 The paper contains charts in "eps" and "png" formats. The numeration is such that 
the "png" charts follow the "eps" ones. 

A quantitative analysis of these properties is presented in the companion paper [5]. 



6 



Thus a simple interchange of price increments Ci <-> C2 is not a symmetry of 
price dynamics. 

The third crucial property of the V(x,y) is that all even sectors are 
"stronger" (containing more probability) than the odd ones. This leads to a 
number of remarkable consequences described in the next paragraph. 

2.2.1 Asymmetry of conditional response 

Let us first analyze the symmetry of the distribution V(x, y) with respect to 
the axis y = 0. In term of sectors shown in Fig. 1 we compare sector I versus 
sector VIII, sector II versus sector VII, etc. To arrange a convenient frame- 
work for analyzing this question let us separate the two-dimensional distribu- 
tion V(x, y) into a sum of symmetric Vf 6 and antisymmetric Vf components: 

V(x,y) = Vf + Pf (2) 

where 

V! = \ (V(x, y) + V(x, -y)) and V? = \ (V(x, y) - V(x, -y)) (3) 

This decomposition is of course not a sum of two probability densities. To 
provide a probabilistic meaning for the asymmetric term Vf(y, x) one should 
consider only its positive part. 

To get a better visual picture of the asymmetric component Vf(x,y) 
let us consider its positive part V^ p \x,y) = 6[Pf(x,y)} Vf(x,y), where 
#[•] is a Heaviside step function. The definition of V a I {:P \x,y) guarantees 
that no information is lost when imposing this restriction. In Fig. 8 we 
show a two-dimensional projection of the surface log 4 Vj ^ (x, y) in the range 
$ — 0.3 <= x, y <= $ 0.3 for the case of two consequent 1 - minute intervals. 

We see that the plot in Fig. 8 demonstrates a remarkable four - blade 
mill - like pattern or simply the market mill pattern. Before analyzing this 
asymmetry in more details, let us point out that the overall asymmetry, if 
measured with respect to the symmetric component, is of order of 5 percent, 
so we are in fact dealing with a small asymmetric component on top of the 
dominating symmetric one. 

To study the structure of the response asymmetry in more details, 

let us consider the sections of its surface by the constant push planes x = 

6 We use a subscript "I" to index the first of the four asymmetries considered in the 
paper. 
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const. In Fig. 9 we show three response profiles for \og 4 Vj for the push 
values of x = $ 0.01, 0.07, 0.25 . 

For small pushes we observe a clear example of the skewed distribution. 
The small negative response after the small initial positive push happens 
quite often, but is outweighted by the positive tail corresponding to big pos- 
itive response. So, a reaction to the small push is either the small contrarian 
response or the big followup in the direction of the initial push separated by 
the gap. At medium pushes the situation becomes more balanced: contrar- 
ian mechanism eats away a large chunk of the response asymmetry space. 
Finally, for large pushes the response is purely contrarian. 

A dependence of the response pattern on the magnitude of the push 
is best quantified by the conditional distribution V(y\x) = V(x,y)/V(x). 
In particular the market mill structure shown in Fig. 8 leads to a very 
nontrivial z-shaped nonlinear behavior of the conditional mean response 
(y)x = / dyyV(y\x). To illustrate this phenomenon we plot in Fig. 2 the 
dependence of the conditional mean response on the push at three time scales 
AT = 1 min., 3 min. and 6 min. 

To clarify the relation between the market mill pattern shown in Fig. 8 and 
the behavior of the conditional response let us note that the decomposition 
of the two-dimensional distribution in Eqs. (|2I3|) generates a corresponding 
decomposition of the conditional distribution V(y \ x) into symmetric and 
antisymmetric parts 

V(y\x) = V 8 I (y\x) + V?(y\x) (4) 

It is easy to check that the conditional response (y) x is a functional of the 
asymmetric component only: 

(y) x = J dyyV{y\x) = J dyyVf{y\x) (5) 

so that the nontrivial nonlinear response shown in Fig. 2 directly follows from 
the market mill pattern of Fig. 8. In terms of sectors shown in Fig. 1 the 
nonlinear mean response follows from the fact that at small x > contribu- 
tion of sector II is stronger than that of sector VIII while at larger x sector 
VIII starts dominating over sector II, etc. 

Let us illustrate the asymmetry of conditional mean by considering a 
sample price trajectory in the xy plane composed by two consecutive price 
increments Ci and £2- The existence of the market mill structure shown in 
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Figure 2: Mean conditional response, 1 min ==>- 1 min, 3 min ==>- 3 min, 6 
min ==>■ 6 min 

Fig. 8 means that the probability ^(Ci, C2) of the combination (Ci ; C2) is not 
equal to that of the combination (d, —(2), see the diagram in Fig. 3. 

How does the market mill asymmetry depend on the length of the time 
interval? To study this question we have examined a set of intervals of 
length AT = 6 and 30 minutes. The corresponding asymmetries are shown 
in Figures. 11-12. All the plots show the presence of the market mill structure 
which is progressively becoming more noisy with the growing length of the 
interval. Let us also mention that if we have a gap between the two time 
intervals under consideration, the market mill structure is still present. This 
is illustrated in Fig. 13 where we show the asymmetry for two 1 minute 
intervals separated by 1 minute gap. 
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Figure 3: Diagram illustrating the conditional response asymmetry 

2.2.2 Asymmetry with respect to increment's ordering 

Let us now analyze the symmetry properties of V(x, y) with respect to the 
axis y = x. The appropriate decomposition of V(x,y) into symmetric and 
antisymmetric contributions reads 

P(x,y)=V II + V! I (6) 

where 

Vu = \ (V(x, y) + V{y, x)) and Pfj = ~ (V(x, y) - V(y, x)) (7) 

A two-dimensional projection of \og 4 V jj P \ where *pff^ is a positive compo- 
nent of the asymmetric part Vfj, is plotted in Fig. 14 showing the market mill 
pattern similar to the one for the asymmetry of the conditional response in 
paragraph 2.2.1. The difference is that for the market mill in paragraph 2.2.1 
the first blade (sector 2) with respect to the axis y = is one sector away 
from it, while here the first blade (the same sector 2) is the sector closest to 
the axis y = x. This difference results in specific z-shaped pattern discussed 
below, see Fig. 4. 

Another way of analyzing the symmetry properties of V(x, y) is to intro- 
duce new coordinates 

1 1 
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which means that we have rotated the original coordinate axes at 4> — 7r/4. 
This specific rotation is chosen because geometrically the symmetry with 
respect to the transformation x < — ► y is the symmetry upon reflection with 
respect to the axis x = y, which in the new coordinates (z, z) a symmetry 
with respect to the axis z — 0. In new coordinates 

P s u = 1 (V(z, z) + V(z, -z)) and V a u = \ (P(z, z) - V(z, -z)) (9) 

The properties of the distribution Q can be conveniently analyzed by con- 
sidering a conditional distribution V(z\ z). In Fig. 4 we plot the conditional 
mean (z) z versus z at 1 min. time scale. We see the z-shaped structure 
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Figure 4: Conditional mean with respect to the axis y = x, 1 min =^> 1 min 

analogous to the one discussed in the previous paragraph, with the differ- 
ence that the conditional mean z is now trend- following at large \z\ > Zq, 
where Zq ~ $0.1, and contrarian at small \z\ < z . This asymmetry of the 
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conditional mean (z) z is, of course, another manifestation of the basic asym- 
metry with respect to the transformation x < — > y. To illustrate conclusions 
following from the pattern plotted in Fig. 4, let us again consider the sample 
path (Ci, C2) and list the basic patterns: 

Ci + C2 < —Zo =>• on average C2 < Ci 

—z < Ci + C2 < =>• on average C2 > Ci 

< Ci + C2 < z Q => on average ( 2 < Ci 

Ci + C2 > zq =>- on average C2 > Ci (10) 

A clear interpretation of the asymmetry with respect to the axis y — x is 
obtained by comparing the probabilities of two sample paths 

• Path 1: Push (1, response (2, Ci + C2 fixed, probability "P(Ci,C2) 

• Path 2: Push ( 2 , response £i, ( 2 + Ci fixed, probability P(C2,Ci) 
This is illustrated in the diagram in Fig. 5. The market mill pattern corre- 




Figure 5: Diagram illustrating the asymmetry with respect to increment's 
ordering 

sponding to the asymmetry of V(x, y) with respect to the axis y = x leads 
therefore to the striking conclusion that the ordering of the price increments 
is probabilistically essential! 
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2.2.3 Conditional push asymmetry 

Let us turn to the analysis of the symmetry of the distribution V(x,y) with 
respect to the axis x — 0. The appropriate separation of the full distribution 
into symmetric and antisymmetric parts reads 

V(x,y) = Vhr + Vfu, (11) 

where 

nu = \(V(x,y)+V(-x,y)) and V a lu = \ (V(x, y) - V(-x, y)) (12) 

We are thus dealing with an "inverse" relation between the push and re- 
sponse, namely, with the structure of push patterns leading to a given re- 
sponse. A two-dimensional projection of log^^ , where vffj is a positive 
component of the asymmetric part Vf n , is plotted in Fig. 15. It shows the 
same market mill pattern as characterizing the asymmetry of conditional 
response discussed in paragraph 2.2.1. All the discussion on the asymme- 
try of conditional response can be repeated here for the conditional push. 
This includes, in particular, the zigzag structure of the dependence of mean 
conditional push on response (x) y . 

Another way of looking at this problem is a rotation of the axes at <fi = 
tt/2: 

z = y ; z — —x (13) 

followed by the corresponding decomposition of the full probability distribu- 
tion with respect to the axis z = 0. 

The asymmetry of conditional push corresponds to considering the prob- 
abilities of the two sample paths (Ci; C2) and (— Ci,C2), respectively P(Ci,C2) 
and V(— Ci, C2), so that we are interested in the properties of the "prehistory" 
of the increment in the second interval ( 2 , see the diagram in Fig. 6. 

2.2.4 Asymmetry of mirror orderings 

Let us finally analyze the symmetry properties of V(x, y) with respect to the 
axis y = —x. The corresponding symmetric and antisymmetric components 
read 

V(x,y) = V s IV + Vf v (14) 
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Figure 6: Diagram illustrating the conditional push asymmetry 



where 

V\ v = ^(V(x,y)+V(-y,-x)) and V a IV = X - (V(x, y) - V(-y, -x)) 

(15) 

A two-dimensional projection of log 4 Vjy \ where Vjy is a the positive com- 
ponent of the asymemtric part Vf v , is plotted for the case of two consecutive 
1 - minute intervals in Fig. 16. We see the same four - blade market mill 
pattern as in paragraph 2.2.2. 

A useful view on the symmetry with respect to the axis y = —x is provided 
by using the rotated coordinate frame, which in this case is obtained by 
rotation at <p = 37r/4: 

z = -^=(x-y); z=-^=(x + y) (16) 

The conditional mean (z) z has the same shape as in the case of the asymmetry 
with respect to the increment's ordering, but now this is a conditional mean 
of the total increment z\ + Z2 at given difference z\ — Z2- 

To elucidate the meaning of the asymmetry with respect to the axis y = 
—x, let us consider the probabilities of two sample paths 

• Path 1: Push £i, response (2, Ci + C2 = z*, probability V((\, C2) 

• Path 2: Push — £ 2 , response — (\, — ( 2 — Ci = ~ z *, probability "P(— C2, — Ci) 
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Figure 7: Diagram illustrating the asymmetry of mirror orderings 

This is illustrated by the diagram in Fig. 7. Note that the second path is a 
mirror one with respect to the first. One could also describe the Path 2 as a 
trajectory reversal. The meaning of the asymmetry with respect to the axis 
y = —x is thus that the probabilities associated with two mirror trajectories 
are not equal. 

2.3 Discussion 

Let us add a few more points to the discussion of the properties of the 
structure of probability distribution in the push - response plane. 

First, the existence of the z-shaped conditional mean response implies 
the existence of a nonlinear predictability 7 . Does this predictability reflect 
market inefficiency at the timescales under consideration? The generic prob- 
abilistic test of the market efficiency is model-dependent fHl 1201 El 1221 123] ■ 
A simple, although not totally generic picture of the efficient market is that 
of a martingale dynamics |2*lj . In this picture the full implementation of 
available information into current stock price means that the expected value 
of this price in the future equals its current one, which in the framework of 
this paper means that (y) x = for all x. The nonzero conditional mean is 
thus explicitly breaking the martingale condition (y) x = 0, see Eq. ©, and 

7 By nonlinear predictability we mean the nonlinear dependence of conditional mean 
response on the push. 
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thus, within the generally accepted paradigm, signals market inefficiency. 

Second, it is of substantial interest to view our results on the condi- 
tional response in connection with an usual measure of dependence between 
stochastic variables - the covariance. The unconditional mean increments 
in all intervals considered are so small that one can safely neglect it. The 
covariance is than equal to 



It is clear the nonlinear structure of the conditional response (y) x strongly 
affects the value of conditional mean. For example, if the conditional re- 
sponse (y) x is even in the push, say (y) x oc x 2 , the covariance and thus the 
linear correlation coefficient are, for symmetric marginal distribution P(x), 
identically zero. A somewhat more realistic example of (y) x = ax — bx 3 gives 
(xy) = a(x 2 ) — b(x i ), so that in this case the value of the linear correlation 
coefficient depends, through parametrization of the conditional response, on 
higher moments of the marginal distributions P(x). We see that generi- 
cally the linear correlation does not provide an adequate description of the 
response. In particular, the zigzag form of the conditional response makes 
description in terms of linear correlation unreliable. Let us notice that even 
the sign of the high frequency autocorrelation does not seem to be well estab- 
lished. The popular view is that at short timescales the increments (returns) 
are anticorrelated j3]. A detailed analysis of high frequency returns in [25] 
showed, however, a positive autocorrelation. This ambiguity can well be re- 
lated to the nonlinear structure of the conditional response discussed in the 
present paper. 

Third, the existence of market mill structure means that price evolution 
is characterized by an explicit breakdown of the symmetry with respect to 
time reversal, which in the considered case is a symmetry with respect to 
the push-response interchange x <-> y. Previously the time reversal asym- 
metry effects were found only in the context of volatility dynamics [26] I2Z] • 
Let us note, that the breakdown of the time reversal symmetry and the 
zigzag structure of the mean response provides valuable information on the 
two-dimensional probability distribution V(x,y). For the popular elliptical 
distributions, i.e. those dependent on the quadratic form dependent on x 
and y, the mean response is always linear in the push and the distribution 
function is symmetric with respect to x *->■ y transformation. Therefore, 
the two-dimensional distribution of consecutive increments is, in fact, non- 




(17) 
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elliptical, which is a serious obstacle in developing an analytical description 
of this distribution. Let us mention that in [28J it was shown that the bi- 
variate distributions of the returns at FX market is not elliptical at short 
intraday timescales. 

The above-described remarkable market structures characterize the bi- 
variate distribution for the whole ensemble of price increments of all the 
stocks considered. Does this blend also give an adequate description for 
asymmetry characterizing individual stocks? Careful analysis shows that in- 
dividual asymmetry patterns vary strongly from security to security. This 
makes the stability in these individual patterns even more unexpected. One 
could say that each stock is characterized by its own "fingerprint" of the 
asymmetry of push-response distribution. A detailed analysis of this issue 
will appear in the forthcoming work [7j. Here we just show in Fig. 17 three 
individual characteristic patterns corresponding to the cases with the domi- 
nance of nonlinear response, ordinary correlation and anticorrelation patterns 
respectively. To illustrate the stability of the asymmetry patterns in Fig. 17 
they are shown in two non-intersecting time intervals. 

3 Conclusions and outlook 

Let us formulate the main results obtained in the present paper. Study- 
ing the full probabilistic description of the dynamics of price increments in 
two nonintersecting intraday time intervals for a set of stocks we found an 
approximate symmetry of the bivariate distribution V(x,y) with respect to 
rotations at the multiples of 7r/2 and four remarkable asymmetries with re- 
spect to the axes y = 0, y = x, x = and y = —x characterized by the 
market mill patterns: 

• The asymmetry of the conditional response. 

• The asymmetry with respect to the increment's ordering. 

• The asymmetry of the conditional push. 

• The asymmetry of the mirror orderings. 

All of four market mill patterns follow from the specific structure of the 
general probability distribution V(x,y). The present study was focused on 



17 



describing and interpreting the market mill patterns having left many ques- 
tions for future analysis jHllZj- Some particular issues we are going to address 
are: 

• A detailed analysis of the shape of response pattern [Hj. Besides giving 
a detailed view on the structure of the response profile this analysis 
provides a direct link to a large body of econometric literature devoted 
to modeling of price dynamics in terms of regressive models such as, 
e.g., that in [H]. 

• We have briefly mentioned the remarkable stability of the asymmetry 
patterns of the individual stocks. A detailed study of these patterns and 
an analysis of the way the noisy individual patterns combine to produce 
a clear overall market mill structure is perhaps the most important 
problem to address [7j. 

• At the fundamental level the main challenge is to build a parsimonious 
model of agent strategies allowing to reproduce the specific structure 
of response asymmetry described in the present paper. In particular, 
it is important to discuss the issue of predictability using a more mi- 
croscopic type of description such as used in jSOj- Work on these issues 
is currently in progress. 

The work of A.L. was partially supported by the Scientific school support 
grant 1936.2003.02 

4 Appendix 

Below we list stocks studied in the paper: 

A, AA, ABS, ABT, ADI, ADM, AIG, ALTR, AMGN, AMD, AOC, APA, 
APOL, AV, AVP, AXP, BA, BBBY, BBY, BHI, BUB, BJS, BK, BLS, BR, 
BSX, CA, CAH, CAT, CC, CCL, CCU, CIT, CL, COP, CTXS, CVS, CZN, 
DG, DE, EDS, EK, EOP, EXC, FCX, FD, FDX, FE, FISV, FITB, FRE, 
GENZ, GIS, HDI, HIG, HMA, HOT, HUM, JBL, JWN, INTU, KG, KMB, 
KMG, LH, LPX, LXK, MAT, MAS, MEL, MHS, MMM, MO, MVT, MX, 
MYG, NI, NKE, NTRS, PBG, PCAR, PFG, PGN, PNC, PX, RHI, ROK, 
SOV, SPG, STI, SUN, T, TE, TMO, TRB, TSG, UNP, UST, WHR, WY 
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Figure captions 



Figure 10. The full bivariate distribution \og 8 V(x, y), 1 min ==>- 1 min. 

Figure 11. The asymmetric component \og 4 V^ p \x,y) with respect to 
the axis y = 0, 6 min ==>- 6 min. 

Figure 12. The asymmetric component log 4 vf p \x, y) with respect to 
the axis y = 0, 30 min ==>- 30 min. 

Figure 13. The asymmetric component log 4 Vj^\x, y) with respect to 
the axis y = x, 1 min ==>- 1 min. 

Figure 14. The asymmetric component log 4 Vjjj(x, y) with respect to 
the axis x — 0, 1 min 1 min. 

Figure 15. The asymmetric component log 4 Vjy \x, y) with respect to 
the axis y = —x, 1 min 1 min. 

Figure 16. The asymmetric component \og 4 Vj^ p \x , y) with respect to 
the axis y — 0, 1 min 1 min, gap between the intervals 1 min. 

Figure 17. The asymmetric components log 4 Vf p \x,y) with respect to 
the axis y — 0, 1 min 1 min, individual stocks. 
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Asymmetry with respect to y=0 axis, 1 min => 1 min 
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Figure 8: Two-dimensional projection of the asymmetry of conditional dis- 
tribution, 1 min ==>■ 1 min 
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Figure 9: Asymmetry of conditional distribution, 1 min 1 min, selected 
push values 0.01,0.07,0.25. 
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